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Abstract 



, The influence of the initial composition of the exploding white dwarf on the nucleosynthesis, light 

■ curves and spectra of Type la supernovae has been studied in order to evaluate the size of evolutionary 
effects on cosmological time scales, how the effects can be recognized and how one may be able to 

^ ■ correct for them. 

CO , The calculations are based on a set of delayed detonation models which give a good account of 

the optical and infrared light curves and of the spectral evolution. The explosions and light curves 
are calculated using a one-dimensional Lagrangian radiation-hydro code including a nuclear network. 
' Spectra are computed for various epochs using the structure resulting from the light curve code. Our 

, NLTE code solves the relativistic radiation transport equations in the comoving frame consistently 

0^ ' with the statistical equations and ionization due to 7 radiation for the most important elements (C, O, 

r-| ! Ne, Na, Mg, Si, S, Ca, Fe, Co, Ni). About 1,000,000 additional lines are included assuming LTE-level 

Qh' populations and an equivalent-two-level approach for the source functions. 

Changing the initial metallicity Z from Population I to II alters the isotopic composition of the 
' outer layers of the ejecta that have undergone explosive O burning. Especially important is the 

c/3 . increase of the ^'^Fe production with metallicity. The influence on the resulting rest frame visual 

' and blue light curves is found to be small. Detailed analysis of spectral evolution should permit a 

determination of the progenitor metallicity. 

Mixing ^^Ni into the outer layers during the explosion can produce effects similar to an increased 

■ initial metallicity. Mixing can be distinguished from metallicity effects by means of the strong cobalt 
' and nickel lines, by a change of the calcium lines in the optical and IR spectra and, in principle, by 

7-ray observations. 

As the C/0 ratio of the WD is decreased, the explosion energy and the '^'^Ni production are 
reduced and the Si- rich layers are more confined in velocity space. A reduction of C/0 by about 60 
% gives slower rise times by about 3 days, an increased luminosity at maximum light, a somewhat 
faster post-maximum decline and a larger ratio between maximum light and "^^Ni tail. A reduction 
of the C/0 ratio has a similar effect on the colors, light curve shapes and element distribution as a 
reduction in the deflagration to detonation transition density but, for the same light curve shape, the 
absolute brightness is larger for smaller C/0. An independent determination of the initial C/0 ratio 
and the transition density is possible for local SN if detailed analyses of both the spectra and light 
curves are performed simultaneously. 

Because the spectra are shifted into different color bands at different redshifts, the effect of 
metallicity Z on a given observed color is a strong function of redshift. A change of Z by a factor of 3 
or the C/0 ratio by 33 % alters the peak magnitudes in the optical wavelength range by up to « 0.3"^ 
for z > 0.2. These variations are comparable to the effect of changes of ftj^.f and A at redshifts of 0.5 
to 1. The systematic effects due to changes in composition are expected to remain small up to about 
z « 0.5 for R-V and up to z « 0.7 for R-I. 
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With proper account of evolutionary corrections, supcrnovac will provide a valuable tool to de- 
termine the cosmological parameters and they will provide new insight into the chemical evolution of 
the universe. 

Subject headings: Supernovae and supernovae remnants: general - hydrodynamics - nucleosynthesis 
- radiation transfer - light curves - Hq, A 



1. Introduction 

Type la Supernovae (SNe la) may reach the same brightness as the entire parent galaxy. In 
principle, this allows for the measurements of extragalactic distances and cosmological parameters 
such as the Hubble constant Hq, f^JVf , fiyV' ^^'^ related deceleration parameter qo- It is now 
widely accepted that Type la supernovae are not a strictly homogeneoiis class of objects with similar 
brightness (e.g. Barbon, Ciatti & Rosino, 1978; Pskovskii 1970, 1977; Branch 1981: Phillips et al. 
1987, Filippenko et al. 1992a; Phillips et al. 1992; Fihppenko et al. 1992ab; Leibundgut et al. 1993; 
Hamuy et al. 1996). For nearby supernovae {z < 0.1), different schemes have been developed and 
tested to cope with the problems of deducing the intrinsic brightness based on theoretical models 
or observed correlations between spectra or light curves and the absolute brightness using primary 
distance indicators (e.g. Norgaard-Nielsen et al. 1989; Branch & Tammann 1992, , Sandage & 
Tammann 1993, Miiller & Hoflich 1994, Hamuy et al. 1996; Riess, Kirshner & Press 1995; Nugent et 
al. 1995; Hoflich & Khokhlov 1996). These methods have been tested locally and provide consistent 
results. New telescopes and observational programs put SNe la at large red shifts well within reach 
and justify optimism for the discovery of a large number of distant SNe la. The Berkeley group has 
discovered more than 50 SNe la up to a red shift of 0.9 (Pennypackcr et al. 1991, Perlmutter et al. 
1997; Pennypacker, private communication). The CfA/CTIO/ESO/MSSSO collaboration has found 
a similar number of supernovae (Leibundgut et al. 1995, Schmidt et al. 1996). Systematic errors due 
to evolutionary effects represent a concern for the use of SNe la to determine the shape of the universe 
(Hoflich et al. 1997). This is especially true for statistical methods which are calibrated only on local 
SNe la. This paper represents a first attempt to characterize some of these evolutionary effects and 
points the way to eliminate associated systematic errors. 

It is widely accepted that SNe la are thermonuclear explosions of carbon-oxygen white dwarfs 
(Hoyle & Fowler 1960; for discussions of various theoretical aspects see Woosley & Weaver 1986, 1995, 
Wheeler & Harkness 1990, Canal 1994, Nomoto et al. 1995, Nomoto 1995, Wheeler et al. 1995, and 
Hoflich & Khokhlov 1996). Three main scenarios can be distinguished: 

A primary scenario consists of massive carbon-oxygen white dwarfs (WDs) with a mass close to 
the Chandrasekhar mass which accrete through Roche-lobe overflow from an evolved companion star 
(Nomoto & Sugimoto 1977; Nomoto 1982). In these accretion models, the explosion is triggered by 
compressional heating. From the theoretical standpoint, the key questions are how the flame ignites 
and propagates through the white dwarf. Several models of SNe la within this general scenario have 
been proposed in the past, including detonations (Arnett 1969; Hansen & Wheeler 1969), deflagrations 
(Ivanova, Imshennik & Chechetkin 1974; Nomoto, Sugimoto & Neo 1976) and delayed detonations, 
which assume that the flame starts as a deflagration and turns into a detonation later on (Khokhlov 
1991, Yamaoka et al. 1992, Woosley & Weaver 1995). The latter scenario, the so-called "delayed 
detonation" and its variation "pulsating delayed detonation" , seems to be the most promising one, 
because, from the general properties and the individual light curves and spectra, it can account for 
the majority of SNe la events (e.g. Hoflich & Khokhlov 1996, Hoflich et al. 1997, Nomoto et al. 1997, 
Nugent et al. 1997, and references therein). In addition, with the discovery of the supersoft X-ray 
sources, potential progenitors have been found (van den Heuvel et al. 1992; Rappaport et al. 1994ab; 
DiStefano et al. 1997). We note that the classical "deflagration" model W7 has a similar structure to 
the "delayed detonation" models which have been successfully applied to reproduce light curves and 
spectra of normal bright SNe la (Harkness, 1991). 

The second scenario for progenitor models consists of two white dwarfs in a close orbit which 
decays due to the emission of gravitational radiation and this, eventually, leads to the merging of the 
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two WDs. In an intcrnicjcliatc step, those models form a low density WD surrounded by a CO envelope 
(Webbink 1984; Iben & Tutukov 1984; Paczynski 1985, Iben 1997). 

A third class of models involves double detonation of a C/O-WD triggered by detonation of 
a helium layer in low-mass white dwarfs as explored by Nomoto (1980), Woosley, Weaver & Taam 
(1980), and most recently by Woosley & Weaver (1994), Livne & Arnett (1995) and Hoflich & Khokhlov 
(1996). From light curves and spectra, this scenario can be excluded as accounting for the majority 
of SNe la events (Hoflich et al. 1997, Nugent et al. 1997). 

A separate but closely related uncertainty is the evolution of the progenitor systems (e.g. Canal 
et al. 1995, Hernanz et al. 1997). If the progenitor population undergoes a change with time, methods 
to determine ^]\,f and A (and hence qo) using only local calibrators could be systematically flawed. 

Time evolution is expected to produce the following main effects: (a) a lower metallicity will 
decrease the time scale for stellar evolution of individual stars by about 20 % from Pop I to Pop II 
stars (Schaller et al. 1992) and, consequently, the progenitor population which contributes to the 
SNe la rate at any given time. The stellar radius also shrinks. This will influence the statistics of 
systems with mass overflow; (b) Evolutionary effects of the stellar population will change the mass 
fmiction present at the time corresponding to a given redshift; (c) The initial metallicity will determine 
the electron to nucleon fraction of the outer layers and hence affects the products of nuclear burning; 
(d) Systems with a shorter life time may dominate early on and, consequently, the typical C/0 ratio 
of the central region of the WD will be reduced; (e) The properties of the interstellar medium may 
change; (f) In principle, a change of the metallicity may alter the M-R relation. Fig. 1 shows that 
this issue is not a major concern. 
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Figure 1: Influence of the metallicity on the structure of a WD. The shaded band marks the region 
which allows for a successful reproduction of observed LCs and spectra for SNe la (Hoflich & Khokhlov 
1996). In the upper plots, the M-R relation and the pc-M relations are given for solar metallicities 
and, in the lower plots, the difference (in %) for models with 1/100 of solar metallicity are given. On 
the scales of the upper plots, the curves corresponding to different metallicities would merge. 
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In this paper wc address the possible influenee of cvohitionary effects on hgfit curves and spectra 
for the delayed detonation scenario. Because the time evolution of the composition is not well known 
and because we have not considered the entire variety of possible models, the values given below do 
not provide a basis for quantitative corrections of existing observations. The goal is to get a first 
order estimate of the size of the systematic effects and to demonstrate how evolutionary effects can be 
recognized in a real data sample and how one may be able to correct for them in the determination 
of cosmological parameters. In §2, the numerical treatment is briefly described. In §3, the influence 
of metallicity and the variation of the C/0 ratio on the isotopic abundances, density and velocity 
structure of the envelope are discussed. Light curves and maximum light spectra are presented. 
Wc discuss briefly how the effects of metallicity can be distinguished from mixing processes. In §4, 
the influence of evolutionary effects on the use of SNe la to determine cosmological parameters is 
presented. Section 5 gives a discussion of our results and conclusions. 

2. Brief Description of the Numerical Methods 

2. 1 . Hydrodynamics 

The explosions are calculated using a one-dimensional radiation-hydro code, including nuclear 
networks (Hoflich & Khokhlov 1996). This code solves the hydrodynamical equations explicitly by 
th{^ pi(X"ewisc! parabolic; mcithod (Collela & Woodward 1984) and includes the solution of the; radiation 
transport implicitly via moment equations, expansion opacities, and a detailed equation of state 
(Hoflich et al. 1993). About five hundred depth points are used. Radiation transport has been included 
to provide a smoother transition from the hydrodynamical explosion to the phase of free expansion. 
We omit 7-ray transport during the hydrodynamical phase because of the high optical depth of the 
inner layers. Nuclear burning is taken into account using our network which has been tested in many 
explosive environments (see e.g. Thielemann, Nomoto & Hashimoto 1996 and references therein). 
During the hydrodynamical calculations, an a-network of 14 isotopes is included to describe the 
energy release. The final chemical structure is calculated by post-processing the hydrodynamical 
model using a network of 216 nuclei. The accuracy of the energy release in the reduced network has 
been found to be about 1 to 3 % . 

2.2. Light Curves 

Based on the explosion models, the subsequent expansion and bolometric as well as monochro- 
matic light curves are calculated using a scheme recently developed, tested and widely applied to SN 
la (e.g. Hoflich et al. 1993, Hoflich et al. 1996 and references therein). The code used in this phase 
is similar to that described above, but nuclear burning is neglected and 7 ray transport is included 
via a Monte Carlo scheme (Hoflich, Miiller & Khokhlov 1992). In order to allow for a more consistent 
treatment of scattering, wc solve both the (two lowest) time-dependent radiation moment equations 
for the radiation energy and the radiation flux, and a total energy equation. At each time step, we then 
use T(r) to determine the Eddington factors and mean opacities by solving the frequency-dependent 
radiation transport equation in the comoving frame (see next section) in about 100 frequency bands 
(see below) and integrate to obtain the frequency- averaged Eddington factors. We use the latter to 
iterate the solution with the frequency-integrated energy and flux equations. The frequency averaged 
opacities, i.e. Planck, Rosseland and Flux means, have been calculated under the assumption of local 
thermodynamical equilibrium. This is a reasonable approximation for the light curve since diffusion 
time scales are always governed by layers of large optical depths. Note that the comparison of L(r) 
between the frequency independent solution and the frequency dependent solution provides a critical 
test for the consistency of the approximations used in the radiation hydro. 

Both the monochromatic and mean opacities are calculated using the Sobolev approximation 
(Sobolev 1957) to calculate the absorption probability within a shell and to include line blanketing. 
The approach is similar to Karp et al. (1977) but generalized for the comoving frame and the integra- 
tion boundaries are adjusted to a radial grid (Hoflich 1990). The scattering, photon redistril)ution and 
thermalization terms used in the light curve opacity calculation are calibrated with NLTE calculations 
using the formalism of the equivalent- two- level approach (Hoflich 1995). 
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Figure 2: Test for the consistency of the 
integrated luminosity based on the so- 
lution of the monochromatic and fre- 
quency integrated radiation transport 
equation for model DD21c (Table 1). 
About 780 frequencies have been used 
in the former case. Note that this in- 
consistency only enters the Eddington 
factors which are based on the solution 
of monochromatic transport equation, 
but not the energy conservation. 

To calculate the monochromatic light curves, we use the T(r) with the time dependence of 
the structure given by the frequency-integrated solution of the momentum equations to solve the 
frequency-dependent transfer equations in LTE every 0.2 to 0.5 days and to get Li, in the observer's 
frame. The broad band light curves are determine by convolution of Lp with the filter functions. 
We use a few hundred frequency bands with the scheme described in the following subsection. Note 
that a proper treatment of the frequency derivatives in the comoving frame equations requires the 
use of about 5 to 10 times more frequencies than frequency bands. To test the consistency between 
the frequency-dependent and frequency-independent calculations we also integrate the Ly and check 
to see how close the resulting L is to that from the solution of the frequency-integrated moment 
equations. The solutions are consistent within 10 % as shown in Fig. 2 (see also Khokhlov et al. 
1993). Monochromatic colors from our light curve code have been compared to colors calculated by 
our detailed NLTE spectral code (Hoflich 1995). Based on this comparison, the solutions for the Type 
la models are good to a few percent near maximum light and deteriorate to 0.4 mag at about 100 
days when we stop the LC calculation. 

2.3. Spectra 

Finally, detailed NLTE spectra have been constructed based on the LC calculations. Thus, the 
effect of energy stored during previous epochs is properly taken into account. The energetics of the 
model are calculated. Given an explosion model, the evolution of the spectrum is not subject to 
any tuning or free parameters such as the total luminosity. A modified version of our code for Nlte 
Extended ATmospheres (NEAT) is used. Although time dependence in the rate equations can be 
included in this code, it has been omitted because it was found to be negligible. For details see 
Hoflich (1990, 1995), Hoflich et al. (1997) and references therein. 

For the NLTE spectra, the density, chemical proflles and the luminosity as a function of the radial 
distance r are given by the hydrodynamical explosion and the LC calculations, including the Monte- 
Carlo scheme for 7-ray transport. The radiation transport equation (including relativistic terms) is 
solved in the comoving frame according to Mihalas, Kunacz & Hummer (1975, 1976, 1982). Blocking 
by weak lines is included in a 'quasi' continuum approximation, i.e. the frequency derivative terms 
in the radiation transport equation are included in the narrow line limit to calculate the probability 
for photons to pass a radial sub-shell along a given direction /j, (Castor, 1974, Abbot and Lucy 
1985, Hoflich, 1990). The statistical equations are solved consistently with the radiation transport 
equation to determine the non-LTE occupation numbers using both an accelerated lambda iteration 
(sec Olson et al. 1986) and an equivalent-two level approach for transitions from the ground state. 
That provides an efficient way to take the non-thermal fraction of the source function into account 
during the radiation transport and, effectively, accelerates both the convergence rate and the stability 
of the systems (Hoflich 1990, 1995). A comparison of the explicit with the implicit source functions 
provides a sensitive tool to test for convergence of the system of rate and radiation transport equations. 

Excitation by gamma rays is included. Detailed atomic models are used for up to the three 
most abundant ionization stages of several elements, i.e. C, O, Ne, Na, Mg, Si, S, Ca, Fe, Co, Ni, 
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taking into account 20 - 30 levels and 80 - 180 transitions in the main ionization stage. Here, we use 
detailed term-schemes for C II, O II, Ne I, Na I, Mg II, Si II, S II, Ca II, Fe II, Co II and Ni II. 
The corresponding lower and upper ions are represented by the ground states. The energy levels and 
cross sections of bound-bound transitions arc taken from Kurucz (1993ab, 1996). In addition to w 
10,000 lines treated in full NLTE, a total of « 1,000,000 lines out of a list of 31,000,000 (Kurucz 
1993a) are included for the radiation transport. For these lines, we assume LTE population numbers 
for each ion. To calculate the ionization balance, excitation by hard radiation is taken into account. 
LTE-line scattering is not taken as a free parameter (Nugent et al. 1997). Instead, LTE-line scattering 
is included using an equivalent- two-level approach, calibrated by the elements treated in full NLTE. 



3. Results 

3.1. Explosion Models 

The influence of the initial metallicity and mixing on light curves and spectra has been studied 
for the example of a set of delayed detonation models with DD21c being the reference model. In 
Table 1, the quantities given in columns 2 to 7 arc: WD mass; pc central density of the WD 
(in lO^g cm^'^), a ratio of the deflagration velocity and local sound speed; p^j. transition density 
(in 10^5 cm~'^) at which the deflagration is assumed to turn into a detonation; C/0 ratio; 

the metallicity relative to solar by mass; Ej^^^ kinetic energy (in IfP^erg); M^j^ mass of ^^Ni (in 
solar units). The parameters are close to those which reproduce both the spectra and light curves 
reasonably well (Nomoto et al. 1984; Hoflich 1995; Hoflich & Khokhlov 1996). 



Table 1: Basic parameters for the delayed detonation models. 



Model 




Pc 


a 


Ptr 


C/0 






Mm 


DD13c 


1.4 


2.6 


0.03 


3.0 


1/1 


1/1 


1.36 


0.79 


DD21c 


1.4 


2.6 


0.03 


2.7 


1/1 


1/1 


1.32 


0.69 


DD23c 


1.4 


2.6 


0.03 


2.7 


2/3 


1/1 


1.18 


0.59 


DD24c 


1.4 


2.6 


0.03 


2.7 


1/1 


1/3 


1.32 


0.70 


DD25c 


1.4 


2.6 


0.03 


2.7 


1/1 


3/1 


1.32 


0.69 


DD26c 


1.4 


2.6 


0.03 


2.7 


1/1 


1/10 


1.32 


0.73 


DD27c 


1.4 


2.6 


0.03 


2.7 


1/1 


10/1 


1.32 


0.69 



Table 2: Total abundances of model DD21c with solar initial composition. 



He 


C 





Ne 


Na 


Mg 


Al 


1.16E-03 


5.44E-04 


5.03E-02 


1.67E-03 


1.05E-05 


1.02E-02 


6.12E-05 


Si 


P 


S 


CI 


Ar 


K 


Ca 


2.06E-01 


2.30E-05 


1.49E-01 


9.50E-06 


3.56E-02 


7.50E-06 


4.00E-02 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


1.02E-03 


2.20E-02 


1.44E-02 


6.87E-01 


1.28E-02 


2.19E-02 


1.61E-05 



To study the influence of a variation of the C/0 ratio, a model (DD23c) with the same parameters 
as DD21c has been constructed but with C/0=2/3 by abundance. The total abundances of the most 
important elements are given for the reference model DD21c in Table 2. Fig. 3 gives the density and 
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Figure 3: Density and velocity as a function of mass for three delayed detonation models (Table 1). 
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Figure 4: Abundances as a function of the final expansion velocity for the three delayed detonation 
models of Figure 3. Both the initial ^"iVi and the final Fe profiles are shown. 
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Figure 5: Abundances of different isotopes as a function of the expansion velocity for models DD21c, 
DD24c and DD27c that have initial compositions of solar, 1/3 solar, and 1/10 of solar, respectively. 

velocity versus mass for the reference model DD21c, a model with a higher transition density (DD13c) 
and a model with reduced C/0 ratio (DD23c). Fig. 4 gives the composition profiles of the major 
elements for the same three models. The overall density and velocity structures arc insensitive to the 
transition density and to the C/0 ratio. Qualitatively, the final burning products can be understood 
by the relation between hydro dynamical and the individual nuclear time scales. The hydrodynamical 
time scale is given by the total energy release during burning of the progenitor. Since most of the 
energy is released by the explosive burning of carbon and oxygen, the energy production per gram 
depends on the initial chemical composition of the WD, i.e. the C/0 ratio. The nuclear time scales 
arc determined by the peak temperature during burning, which depends on the energy release per 
volume because the energy density is radiation dominated. As the energy release per gram is fixed, 
the peak temperature is given by the composition and the local density. Thus, the latter two are the 
dominant factors for the final composition of a zone. With decreasing transition density, less ^^Ni is 
produced and the intermediate mass elements expand at lower velocities because the later transition 
to detonation allows for a longer pre-expansion of the outer layers (DD21c vs. DDI 3c). Similarly, 
with decreasing C/0 ratio in the progenitor, the specific energy release during the nuclear burning is 
reduced (DD23c vs. DD21c) and the transition density is reached later in time, resulting in a larger 
pre-expansion of the outer layers and a more narrow region dominated by Si. Moreover, a change of 
the C/0 ratio from 1 to 2/3 reduces the ^^Ni production by about 20 %. The kinetic energy is reduced 
by 10 % corresponding to a 5 % change in the mean expansion velocity. To first approximation a 
reduction in the transition density has a similar effect than a decrease of the C/0 ratio or a reduction 
of the central density (Hoflich & Khokhlov 1996). 

To test the influence of the metallicity for Z > 20 (i.e. nuclei beyond Ca) we have constructed 
models with parameters identical to DD21c but with initial metallicities between 0.1 and 10 times 
solar (Table 1). The energy release, the density and velocity structure are virtually identical to that 
of DD21c. The main difference is a slight increase of the 'J"Ni mass with decreasing metallicitv due 
to a higher Ye- The reason is that the metallicity mainly effects the initial CNO abundances of a 
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Tabic 3: Isotopes Mj in Mq with Mj > i.E — lo for delayed detonation models (see Table 1). 



Isotope 




DD21c 




DD13c 




DD23c 




DD24c 




DD25c 




DD26c 




DD27c 




1, 


. IDrj- 


-Uo 


7, 


"^p 


-Uo 


2, 


.Oorj- 


n^ 
-Uo 


1, 


1 fip 

.lOll;- 


n^ 
-Uo 


1, 


1 fip 

.lOll;- 


n^ 
-Uo 


1, 


1 fip 


nQ 
-Uo 


1, 


.16E-03 


12 c 


5, 




D/l 


3 


^np 


-U4 


1, 


n/ip 


n"? 

-Uo 


5, 


.OOIl;- 


n/i 

-U4 


5, 


1QP 


CiA 
-U4 


5, 


/inp 

.4U11;- 


V\A 
-U4 


5, 


.33E-04 


15 N 


6 


.OZH/- 


no 
■uy 


7, 


9np 


■uy 


8, 


nop 


nn 
■uy 


1, 


7SP 


ns 

■Uo 


2, 


1 SP 
. loHj- 


nQ 
■uy 


4, 


fifip 
.DDH/- 


ns 

■Uo 


5, 


.38E-10 


16 Q 


5 


.Uoll/- 


n9 
■uz 


4, 


.DOH/- 


(19 

■uz 


6, 


snp 


n9 
■uz 


4, 


QSP 

.yoU;- 


n9 
■uz 


5, 


ni p 

.UlU;- 


n9 
■uz 


4, 


QQP 

.yyn/- 


n9 
■uz 


4, 


.88E-02 


20Ne 


1, 


.0 / rj- 


-Uo 


1, 


.04rj- 


-Uo 


2, 


npip 

.UOHj- 


n'? 

-Uo 


1, 


fi7P 
.0 ( Ej- 


n"? 

-Uo 


1, 


fif;p 


n"? 

-Uo 


1, 


fiSP 

.OoH;- 


nQ 

-Uo 


1, 


.61E-03 


21 Ne 


1 


.Z i H;- 


-Uo 


6, 


.OOllr 


no 
uy 


2, 


99P 


nc 

-Uo 


1, 


1 fiP 


nc 

-Uo 


1, 


A QP 
.4oJl;- 


ns 

-Uo 


1, 


1 9P 
. IZJl/- 


ns 

-Uo 


1, 


.82E-08 


22Na 


2 


91 TP 


■Uo 


1. 


1 "ip 


■Uo 


4. 


1 SP 


na 

■Uo 


2. 


■^QP 

.oynj- 


ns 

■Uo 


1. 
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star. These are converted during the pre-explosion stellar evolution to N in H-burning and via 
^^N(a,7)l^F(/3+)1^0(a,7)22Ne to nuclei with N=Z+2 in He-burning. The result is that increasing 
metallicity yields i.e. a smaller proton to nucleon ratio Ye throughout the pre-explosive WD (Thiele- 
mann et al. 1997). Higher metallicity and smaller Ye lead to the production of more neutron- rich Fe 
group nuclei and less Ni. For lower metallicity and. thus, higher Ye, some additional ^"Ni IS pro- 
duced at the expense of '^'^Fe and '-'^Ni (Thielemann, Nomoto & Yokoi, 1986). The temperature in 
the inner layers is sufficiently high during the explosion that electron capture determines Ye- In those 
layers, the initial metallicity has no influence on the final burning product. The total production of 
isotopes of all the models is given in Table 3. Fig. 5 gives the distribution of various isotopes in models 
DD21c, DD24c, and DD27c which represent solar, 1/3, and 10 times solar metallicities, respectively. 
The main differences due to changes in Z are in regions with expansion velocities in excess of w 12000 
km/sec. Most remarkable is the change in the ''^Fe production which is the dominant contributor 
to the abundance of iron group elements at these velocities since little cobalt has yet been decayed 
near maximum light. For increasing metallicity, the production of this isotope increases significantly. 
For 1/3 solar metallicity, hardly any Fe is produced at high velocities, but ^^Fe is as high as 5 
% by mass fraction if we start with ten times the solar metallicity. We note that these layers with 
v> 12. 000 km, dominate the spectra around maximum light and lines of iron group elements are 
an important contributor to the line opacities (see below). 

3.2. Light Curves 

As discussed in previous papers (e.g. Hoflich et al. 1997 and references therein), bolometric and 
monochromatic light curves provide a valuable tool to probe the underlying explosion models, namely 
the absolute amount of ^^Ni and its distribution. 



DD21C 
DD23C 
DD24C 



43 

Figure 6: Comparison of bolomet- ^ 
ric light curves of the delayed det- ^ 
onation models DD21c, DD23c and » ''^'^ 
DD24c with otherwise identical pa- 
rameters but with different C/0 ra- 
tios and metallicity relative to so- 
lar (C/0; %) of (1;1), (2/3;l) and 
(1;0.3), respectively. 

time [days] 

Bolometric and monochromatic LCs are shown in Figs. 6 and 7. Reducing the ratio C/0 from 
1/1 to 2/3 in the WD reduces the ^^Ni production and the kinetic energy in DD23c compared to 
DD21c. Both the bolometric and monochromatic light curves are affected. The smaller expansion 
due to the smaller -Ey^j^ causes a reduced geometrical dilution of the matter and a reduction of the 
expansion work at a given time. Consequently, for reduced C/0, the rise time to maximum light is 
slower by about 3 days and the maximum brightness is larger because more of the stored energy goes 
into the radiation rather than kinetic energy. The reduced heating at the later time of maximum 
light also impHes that DD23c is shghtly redder (A(B - V) ^ 0.02"^) than DD21c at its peak. After 
about day 35, however, the luminosity is smaller by about 10 % for lower C/0 due to the smaller 
^^Ni production (i.e. the instantaneous energy production). The postmaximum decline is steeper. 
The change of the C/0 ratio from 1/1 to 2/3 has a similar effect on the colors, light curve shape and 
the distribution of elements as a 10 % reduction of the transition density or of the central density in 
delayed detonations. For a given ratio of peak to tail, a model produced by a reduced C/0 ratio will 
be brighter than a one obtained by varying the transition density (Hoflich 1995). The rise time and 
the expansion rate as measured by the doppler shift of lines (see below) provide a way to determine 
the C/0 ratio and the transition density independently (Table 2, compare Khokhlov et al., 1993, 
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Figure 7: Comparison of light curves in B (left) and V (right) of the delayed detonation models 
DD21c, DD23c and DD24c with otherwise identical parameters but different different C/0 ratios and 
metaUicities relative to solar (C/0; R^) of (1;1), (2/3;l) and (1;0.3). The difference relative to DD21c 
and the monochromatic LCs are given in the upper and lower plots, respectively. 



Hoflich 1995, Hoflich & Khokhlov 1996). The current state of the art for LC calcidations and for the 
progenitor evolution does not allow for a sufficiently fine discrimination to determine the absolute C/0 
ratio. A careful differential analyses of observations may yield information on the relative variation if 
the C/0 ratio, and the central ignition density or transition density. 

Changing the initial mctallicity Z has very little influence on the bolometric light curve (compare 
DD21c, DD24 in Fig. 6) because the ^^Ni production and energy release vary by only 4 % over the 
entire range of models. In addition, diffusion time scales are mainly determined by deeper layers which 
are unaffected by Z. 

Monochromatic LCs arc slightly more affected by changes in the initial metallicity Z (Fig. 7, 
upper panel). The principle effect is that the radiative cooling in the outer layers increases with Z. 
Consequently, depending on the phase, the B-V is larger are redder by about 0.02 to 0.04™. The 
absolute brightness varies by about the same amount. The time of maximum light shifts by only « 1 
day. 

3.3. Spectra at Maximum Light 

Variations in the pattern of the most abundant elements is similar for changes in the C/0 ratio, 
the transition density and the central density of the WD (see Figs. 3 and 4). Consequently, the 
variation of the spectra as a result of changing these parameters is similar. Differences would show 
up in correlation between light curve shape and Doppler shifts of lines at a given phase (sec above). 

The influence of the initial metallicity Z is interesting primarily because it influences the iron 
group elements. As an example, the spectra of the delayed detonation models with solar and 1/3 solar 
metallicity are given in Fig. 8. At maximum light (w 17 days), the line forming region (r « 0.1 ... 1.) 
extends between 1 and 2x10"'^^ cm in the optical (Fig. 8, lower graph), corresponding to expansion 
velocities between 8000 and 16,000 km s^^ . In this velocity range the material has been subjected to 
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Figure 8: Comparison of synthetic NLTE spectra at maximum light for initial compositions of solar 
and 1/3 of solar, respectively (upper graph). The standard Johnson filter functions for UBV, and R 
arc also shown. In the lower graph, the radius as a function of wavelength where the monochromatic 
optical depth reaches 0.1 and 1, respectively is given for DD21c. 

at least partial oxygen burning to silicon and the deeper layers are rich in ^"Ni. This means that the 
iron peak elements visible in the optical spectrum near maximum light are freshly synthesized and do 
not directly represent the initial metallicity. The initial metallicity does affect the optical spectrum 
near maximum, but only indirectly and only in the blue. The abundances of the intermediate mass 
elements are not affected by the metallicity, but the abundance of "^^Fe is a sensitive function of the 
initial Z and "^^Fe provides an important source of Fe lines to the opacity for the matter with v ^ 
10,000 km s-1. The iron line opacity is most substantial for wavelengths less than ~ 4000 A. The 
effect of the initial metallicity on the spectrum can be seen in the wavelength region ^ 3000 A in Fig. 
8 (top). The initial metallicity might be directly observed by examination of UV spectra that form in 
the outermost layers with v ^ 16,000 to 17,000 km s"-*- that have not been subject to oxygen burning. 

By 2 to 3 weeks after maximum, the spectra are completely insensitive to the initial Z because 
the spectrum is formed in even deeper layers where none of the important abundances are affected 
by the metallicity. Thus, for two similar bright SNe with similar expansion velocities, a comparison 
between the spectral evolution can provide a method to determine the metallicity difference or may 
be used to detect evolutionary effects for distant SNe la (see below). 

Note that the UV may provide a unique tool to probe the very outer layers which may have 
undergone little burning during the explosion. 

3.4- Mixing during the Explosion 

Up to now, we have omitted a potential difficulty involved in the determination of Z. As already 

mentioned, we may see strong iron group elements in the line forming region near maximum because 
they are produced at high velocities. Alternatively, these elements may be mixed out during the 
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explosion but after burning took place. To test this, we have artificially mixed the composition but 
kept the overall structure fixed (see Fig. 9). Unlike a change in the initial metallicity, mixing alters 
the entire spectrum. All burning products will be enhanced in the outer region including Ni, Co, and 
Ca. In particular, the Ca lines at about 4000 A and in the IR triplet provide a way to distinguish 
mixing from metallicity effects. 



Fell/CoIISII solar400mix- 




4000 6000 8000 
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Figure 9: Comparison of synthetic spectra at maximum light assuming complete mixing of the outer 
layers, involving 64 % (solar200mix), 48 % (solar300mix) and 24 % (solar400mix) of the total mass, 
respectively (see Figs. 3 & 4). 



4. Implications for Cosmology 

Systematic effects that must be taken into account in the use of SNe la to determine cosmological 
parameters include technical problems, changes of the environment with time, changes in the statistical 
properties of the SNe la, and changes in the physical properties of SNe la which are the main subject 
of this study. The estimates of the redshifts when evolutionary effects become important given below, 
are based on the assumptions (!) of a flat universe with an age of 14 Gyrs, and of evolution time 
scales for the progenitors of about 0.5 to 6 Gyrs which are given by the ZAMS life times of Pop II 
donor stars with masses between 1 and 7 Mq (Schaller et al. 1992). 

In the first class of technical problems, we put corrections for redshift. If standard filter systems 
are used, they can be well calibrated to local standards but the k-correction is of some concern. This 
problem can be overcome if redshifted "standard" filters are used. Another technical problem may 
arise then from the fact that the transmission functions must be identical to those resulting from the 
redshift because no direct calibration can be applied by using a comparison star. 

In the second class of problems associated with enviroment, an important example is that the 
properties of dust may change at high redshift. In the first place, the element abundances in the ISM 
can change. In addition, important donors of dust such as low mass red giant stars cannot contribute 
because their evolutionary time is comparable to or longer than the age of the universe at z w 0.5 to 
1. Another problem related to the correction for extinction is that the extinction law to be applied 
depends on the redshift of the absorbing dust cloud. For details, see Hoflich & Khokhlov (1996). 

In the third class of problems related to the statistical properties is the fact that the contribution 
of different progenitor types may change with redshift. For local supernovae, it is likely that a large 
variety of binary star properties (total mass of donor stars, separation, etc.) with very different 
evolutionary life times account for the variety of SNe la observed (Hoflich et al. 1997). Given this 
variety, we expect a time evolution of the statistical properties of the progenitors because, early on, 
progenitors with a short life time will dominate the sample. Among other changes, more massive, 
shorter lifetime progenitors will have a lower C/0 ratio in the center (see below). Changes in the 
statistical properties are expected to increase substantially at redshifts larger than 0.7 to 0.8 where 
the age of the universe becomes comparable to or shorter than the suspected progenitor life times. 
Other evolutionary effects have been mentioned in the introduction: The ZAMS life time changes 
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with motallicity, Roche lobe radii change with metaUicity and the lower hmit for accretion and steady 
burning of hydrogen on the surface of the WD changes by a factor of 2 between Pop I and Pop II 
stars (Nomoto et al., 1982). 

FinaUy. the properties of the typical, individual SNc la may change. If more massive stars 
contribute to the supernova population, we expect a smaller C/0 ratio. Changing the C/0 ratio from 
1/1 to 2/3 with otherwise identical parameters (see §3.2), will result in a smaller ^^Ni production and, 
consequently, a lower bolometric luminosity at late times. On the other hand, the slowc;r expansion 
causes less adiabatic cooling during early times and, thus, the luminosity at maximum light is larger 
(Fig. 6). This implies that the peak to tail luminosity ratio changes with the C/0 ratio. The 
consequence for analyses that use light curve shapes or the postmaximum decline (e.g. 5mie^) is 
evident. Quantitatively, for monochromatic light curves in our example, the difference in the peak to 
tail luminosity ratio is « 0.1 ... 0.2"^. For techniques to determine the absolute peak luminosity by 
measuring either Smit^ or LC shapes (e.g. Hamuy et al. 1995, Riess et al. 1995), this translates into 
a systematic error of « 0.3"^ in mv. This systematic change in peak luminosity with redshift due 
to changing C/0 ratios could mimic the effects of cosmology. Note that the transition density from 
deflagration to detonation may depend on the energy release during the deflagration phase and hence 
the transition density could be a function of the C/0 ratio. Both effects could alter the LC. 

The other physical effect on the SNe la themselves is the influence of the initial metallicity Z on the 
nuclear burning conditions during the explosion. This produces a change of the isotopic composition in 
the outer layers. A reduction of Z with redshift z is expected. Note that the metallicity corresponds to 
the redshift when the progenitor is formed and not the redshift when the SN la is observed. Although 
small, the changes in the spectrum with metallicity have an important effect on the colors of SN at 
high red shifts where they are shifted into other bands (see Fig. 10). For local SNe la, a change of the 
metallicity by a factor of 3 implies a variation in color of two to three hundredths of a magnitude. For 
z > 0.2, however, Fig. 10 shows that such a change can induce a systematic change in B-V of up to 0.3 
magnitudes. For V-R and R-I, the effect remains small for redshifts z < 0.5 and 0.7, respectively, but it 
increases to the same order at redshifts relevant for the determination of cosmological parameters. The 
amplitude of this effect and hence the uncertainty in color, reddening and brightness with metallicity 
is again comparable to the brightness change imposed by cosmological deceleration. 
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Figure 10: Metallicity effect on B-V and R-V for DD21c as a function of the redshift. 

If SNc at moderate k, 0.5 and high red shifts (« 1.0) are used to determine Vlj^ and O^, 
respectively, these systematic shifts and changes of the amplitude of the variation the influence of the 
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nictallicity on the color indices may prevent a proper determination of the cosmological parameters 
if it is not taken into account. Due to the effect of metaUicity on color and the different time scales 
for progenitor evolution and hence statistical properties discussed above, the distribution of peak 
brightness as a function of light curve shape and redshift is likely to produce a scatter in the distribution 
of the cosmological parameters. The way to reduce this scatter is to understand its physical origin in 
the sorts of effects we have outlined here. Detailed spectral analyses of different SNe may show that the 
effects of metallicity remains moderate (e.g. because, even today, all progenitors may have a low Z), 
or that it can be corrected by using models. Quantitatively, other dynamical models and metallicities 
will give different amplitudes of A(i3 — V), A{V — R), and A{R — I) at given z than illustrated. 
In general, the use of red shifted filters will eliminate the major systematic effect illustrated in Fig. 
10 but, even putting aside the systematic problems due to the realization of identical transmission 
functions, the small 'noise' (see Fig. 10) due to the 'bumpy' spectra is real and requires a very narrow 
grid of filters with respect to redshift, or the determinations of qo will show a significant, intrinsic 
spread. Even if corrections for redshift and the small wiggles are applied and if corrections for the 
interstellar redding are determined, a systematic effect of the order of O.l"^ will remain because of 
the small, but systematic dependencies in the color indices at z=0 (see Fig. 10). 



Table 4: Uncertainties in qo, and i?A as introduced by composition effects 
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In conclusion, both the change of Z and the C/0 ratio in the progenitor may each be expected to 
give variations in the peak brightness derived of the order of a few tenth of a magnitude as one goes 
from z« to z«l. As a guideline, we give in Table 4 the size of the systematic effects as they enter 
qo , and flj^ to first order for Friedmann-Lemaitre cosmological models (Goobar & Perlmutter 

1995, Perlmutter et al. 1997) for which qo = ^ji^ /2 — flj^ with flj^ = A/(3 Hq) and, for a flat univers, 
~^ ~ ^- ^^^^ order, the uncertainties are oc 1/z and Am. An uncertainty Am = 0.6"^ 
should be regarded as worst case scenario. At this point, a warning is appropriate. Although the 
assumed size of the change in the initial metallicity Z is of the right order based on the chemical 
evolution of galaxies (Samland et al. 1997) metallicity change as a function of galactic evolution is 
highly uncertain. Yet another uncertainty enters because Z is set by the time the progenitor is born 
and there can be a large range of times before the star dies. Moreover, the statistical properties of 
progenitor systems are not well known. For instance, the mass of the progenitors and the donor stars 
must be in the range between « 1 and 7 Mq, but the relative contributions within this mass range 
are unknown even for local SNe la. To beat down these uncertainties, comparative analyses of the 
light curves and spectra (e.g. velocities, lineshifts, spectral flux distributions) of both the local and 
distant SNe la are required. 

5. Final Discussion and Conclusions 

Using a delayed detonation model as an example, we have studied the possible influence of the 
initial chemical composition on the light curves and spectra of SNe la. 

At early epochs the mean life time of progenitors will be smaller than at present and, hence, 
the mean progenitor mass is larger compared to the current epoch (see §4). A substantially reduced 
C/0 ratio in the inner region of the progenitor white dwarf will result if the WD originates from a 
star with more than w 3 — 4Mq rather than 1-2 Mq on the ZAMS. We find that a reduction of 
the initial C/0 ratio has effects similar to a lower transition density from deflagration to detonation. 
In both cases, the location of the Si-rich material is shifted to smaller expansion velocity. In the 
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case of a rcidiiction in C/0, the expansion velocity during the explosion is somewhat lower and the 
time is delayed until the transition density is reached. For lower transition densities, the extent of the 
deflagration phase is increased because a longer time is required to reach the transition density. In both 
cases, this allows for a longer pre-expansion of the outer layers and, consequently, the ^^Ni production 
is reduced. The major difference between a reduced C/0 ratio and a lower transition density is the 
reduced explosion energy (and expansion rate) in the former case. This causes the peak in the Si 
distribution to be more pronounced in velocity space, a slower rise to maximum light by about 3 
days, an increased peak luminosity and, consequently, a steeper decline after maximum light when 
the stored, thermal energy is exhausted. The use of light curve shapes or post-maximum decline rates 
alone may lead to systematic errors in the estimates of the absolute brightness by several tenths of a 
magnitude. In principle, detailed coordinated analyses of spectra and light curves provide a way to 
disentangle the effects e.g. of a change in the transition density and the C/0 ratio and to determine 
the peak brightness on a case by case basis. Empirically, a correlation analysis of the LC shape and 
the expansion rate provides a tool to test for the range of C/0 ratios realized in nature (see sect. 3.2). 

Changing the initial metallicity from Population I to II reduces the progenitor life time and Roche 
lobe radii and hence changes progenitor evolution. It also alters the isotopic composition of the outer 
layers. Especially important is the increase in the ^^Fe production with metallicity. The initial WD 
composition has been found to have rather small effects on the overall LCs. The Ni production and 
hence the bolometric and monochromatic optical and IR light curves differ only by a few hundredths 
of a magnitude. This change is almost entirely due to the small change in the Ni production and 
not due to a change in the opacities because the diffusion time scales are governed by the deeper 
layers where burning is complete. The time of maximum light in the Johnson filter system changes 
by w 1 day as Z is varied from solar to 1/3 solar. The short wavelength part of the spectrum at 
maximum light is affected by a change in Z. This provides a direct test for the initial metallicity of 
local SNe if well calibrated spectra are available and, thus, may give a powerful tool to unravel the 
nature (and lifetime) of SNe la progenitors. To do so, spectra around or prior to maximum light 
are needed because at later times the inner layers dominate the spectra where the spectrum is not 
sensitive to the initial metallicity (see sect. 3.3). Note that quantitatively, the influence of Z is model 
dependent because it depends on the photon redistribution and optical depth at a given time, i.e. on 
the density and abundance structure. Even the sign of the effect if the configuration is very different 
than the delayed detonations. As an empirical test, a differential comparison between maximum and 
post-maximum spectra of different SNe la with similar light curve shapes and expansion velocities 
provides a valuable tool to probe for Z effects (see sect. 3.2). 

In terms of distant SNe, we may expect less metals and lower C/0 ratios in the past. For red- 
shifted supernovae, the Z effect may produce systematic changes of up to 0.2 to 0.3 magnitudes in 
B-V, V-R and R-I for z > 0.2, 0.5 and 0.7, respectively. The amplitude of the maximum change and 
the phase with respect to z is somewhat model dependent. In general, the C/0 ratio will decrease with 
redshift and can be expected to become significant at redshifts > 0.7 when the age of the universe is 
comparable to the lifetime of the progenitors. The systematic effects due to changing Z and C/0 may 
be comparable in amplitude to the effects due to the deceleration of the universe. Other evolutionary 
effcicts have been qualitatively discussed in §4 and the introduction. B(X'ause the properties of SNe la 
may change with z, a measurement of qo, ^]<if and A may be highly biased. Determinations of Hp 
will remain unaffected if it is based on a sufficient large range of light curves (e.g. B, V, R, I). 

Finally, we also want to mention the limitations of this study which need to be overcome in the 
future. First of all, only one specific set of models has been studied in detail. Although the model 
parameters have been chosen to allow for a representation of "typical" SNe la, more comprehensive 
studies and detailed fitting of actual observations are needed. The C/0 ratio has been changed for the 
entire WD; however, the C/0 ratio is expected to change only for the mass before binary accretion 
begins (e.g. 0.6 < < I.2M0). AC/0 ratio of « 1 will result from the accreted material at the 
surface regardless of the mass of the primary of secondary star (e.g. Nomoto et al. 1984). The final 
composition of the supernova ejecta which undergo partial burning depends mainly on the density 
during burning. This, in turn, depends on the pre-expansion during the deflagration phase because it 
lasts longer by an order of magnitude (i.e. about 1-2 seconds) compared to the detonation phase. 
During the deflagration phase, only the inner regions (~ 0.2 — 0.3Mq) are burned. Therefore, it does 
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not make imich difference whethcii' the C/0 ratio is lower on the inside; and 1:1 on tlic outside due to 
accretion whether a globally reduced C/0 ratio is considered. Nevertheless, a study of the influence 
of realistic structures is desired in order for the effects to be analyzed more quantitatively. 

This study explored some effects that can be expected if we go back in time. Detailed comparisons 
with observations are ultimately needed to address the effects of composition evolution or to probe 
for their unimportance. To do so, well-sampled spectral and light curve data must be available such 
as currently obtained by the Berkeley and CTIO groups. To isolate the Z effects we would suggest 
that maximum spectra and postmaximum spectra be taken routinely. To acquire a full understanding 
of SNe la, detailed analyses of nearby (e.g. Virgo or closer) supernovae are needed. Such a program 
is justified, because there is already evidence in the current supernova sample for evolutionary effects 
that are not yet understood, but suggest that some effects of varying composition effects are present. 
We might naturally expect a range in initial metallicities of SNe la within spiral galaxies that have 
ongoing star formation and between ellipticals and spirals. The same is true for the C/0 ratio which 
will depend on the average age of the progenitor population. Branch et al. (1996) and Hamuy et 
al. (1996) have shown that the mean peak brightness is dimmer in ellipticals which lack a young 
population. Wang et al. (1997) have shown that the peak brightness in the outer region of spirals is 
similar to those found in ellipticals but, in the central region, both intrinsically brighter and dimmer 
SNe la occur. Both these effects may involve a change in the C/0 ratio. SNe la with a low C/0 ratio 
have a short progenitor lifetime and are predicted to be brighter. Another independent statistical 
test for the influence of a changing C/0 ratio may be to use the spread in the velocity, light curve 
shape relation. If the C/0 ratio is the dominant effect to explain the spread in the properties of LCs, 
there should be a tight relation between LC shape and expansion velocity in ellipticals and a spread 
in spirals. Another approach to get insight into the local SNe la may be the use of J-Cephei distances 
to evaluate the spread in the absolute brightness. 
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